Thermodynamic consistency and fast dynamics in phase field crystal modeling 
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A general formulation is presented to derive the equation of motion and to demonstrate thermo- 
dynamic consistency for several classes of phase field models at once. It applies to models with a 
conserved phase field, describing either uniform or periodic stable states, and containing slow as 
well as fast thermodynamic variables. The approach is based on an entropy functional formalism 
previously developed in the context of phase field models for uniform states [P. Galenko and D. Jou, 
Phys. Rev. E 71, 046125 (2005)] and thus allows to extend several properties of the latter to phase 
field models for periodic states (phase field crystal models). In particular, it allows to demonstrate 
the concept of thermodynamic consistency for phase field crystal models with fast dynamics. 

PACS numbers: 05.70.Fh: 05.70.Ln: 64.60.-i 



Phase field modehng has become a versatile tool for 
studying the dynamics of systems out of equihbrium, 
and has found numerous applications in materials sci- 
ence pj, I2I, [sl, y . For instance, consider a material that is 
disordered at high temperature and has two stable phases 
at low temperature. Upon quenching the material from 
high to low temperature, grains with the different stable 
phases will develop, grow and compete with each other. 
Phase field modeling is able to describe the time evolu- 
tion of such a process. In a phase field model for this 
example, a continuous function of space and time ^'(x, t) 
is introduced - the phase field variable - that assumes 
a different constant value for both stable phases. Near 
an interface between two grains, the value of ^ rapidly 
changes. The phase field variable in this example can be 
interpreted as an order parameter to represent the rela- 
tive mass fraction of both phases. A major advantage of 
phase field modeling is that the phase field contains in- 
formation on the location of all interfaces in the system, 
without the need for explicit interface tracking. Another 
advantage is that it focusses on gross features that are 
common to many systems. System-specific details are in- 
corporated by the interpretation given to the phase field 
variable, and by the values given to the phcnomenolog- 
ical parameters in the model equations. Systems that 
are similar according to appropriate criteria, are there- 
fore described by the same phase field model. Several 
important classes of such models are ordered in Tab. [D 
The major goal of this paper will be to present a unified 
description of four of these classes. 
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TABLE I: A classification of eight different phase field mod- 
els according to the nature of the equilibrium state they de- 
scribe (uniform or periodic) and the properties of the phase 
field variable (conserved or non-conserved) . Models that de- 
scribe slow dynamics only are distinguished from models that 
include fast dynamics as well. For each model, the name of 
the model and the equation of motion is given, as well as the 
free energy (where appropriate). The acronyms refer to Cahn- 
Hilliard (CH), Allen-Cahn (AC), Phase-Field Crystal (PFC) 
and Swift- Hohenberg (SH) - see text. The M in the last line 
stands for 'Modified' throughout. This paper presents a uni- 
fied description for the models that are printed in underlined 
bold face, and demonstrates thermodynamic consistency for 
all four of them. 
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One gross system property to distinguish different 
classes, is whether its stable states are spatially uniform 
or periodic (Tab. |T]). If we neglect details at the atomic 
scale, the two stable states of our illustrative quenched al- 
loy can be considered as being spatially uniform. On the 
other hand, Rayleigh-Benard convection in a thin fluid 
layer between two plates leads to periodic stable states 
with regular patterns of convection cells [l5| . A free en- 
ergy functional Fi[^p\ is postulated for the uniform and 
periodic cases, respectively (Eqs.[l]and[2]). From this free 
energy functional, an equation of motion for the system 
can be obtained by Eqs. [3] or H) Eq. [3] applies when the 



integral of the phase field variable over space is time inde- 
pendent, in which case the phase field is called conserved 
(Tab. IJ). In the other case, Eq. |3]appHes and the phase 
field is non- conserved. The distinction between conserved 
and non-conserved order parameters is another criterion 
to distinguish difi'erent classes. The fact whether the 
equations of motion properly describe the system under 
consideration or not, is an a posierion justification of the 
postulate for the free energy functional. 

Eq. ([T|) was originally introduced by Landau, and since 
then became the fundamental functional to derive phase 
field (PF) models. The potential f{ip) has a double- 
well structure corresponding to the two stable phases de- 
scribed by the order parameter ip. The squared gradient 
term jVV'P accounts for fluctuations in the order param- 
eter and e is a phenomenological parameter that needs 
to be tuned for every particular system. Systems with 
uniform stable states and a conserved order parameter 
were studied half a century ago [5|, and their equation 
of motion is now called the Cahn-Hilliard equation. The 
quenched alloy mentioned before is an example of this. 
Uniform equilibrium states with a non-conserved order 
parameter are described by the AUen-Cahn equation [6|. 
Further examjales have beeii given in the context of poly- 
mer mixtures 
in cosmology 

For systems with stable states that are not spatially 
uniform, the Landau form of the free energy functional 
Eq. ^ is not valid any more. The appropriate energy 
functional, written in a general form IJ], is now Eq. [21 
The symbols have the same meaning as for Eq. [l] only 
there are now two phenomenological parameters ei and 
£2- The Swift-Hohenberg equation [15[ describes systems 
where the stable states are periodic and the order param- 
eter is non-conserved. For the case of periodic systems 
with a conserved order parameter, an equation of motion 
has only recently been derived: the Phase Field Crys- 
tal (PFC) model [l3|, ll4|- It is a simple version of more 
elaborate density functional theories of liquid/crystal in- 



alloys [8|, lof , liquid-crystals 10|, [ill , and 



tcrfaces [16|, |17[. This model captures features at the 



atomic scale, and thus contains highly detailed physical 
information about the system's structure. Such models 
can describe many of the basic properties of polycrys- 
talline materials that are realized during non-equilibrium 
processing. 

An important concept in phase field modeling is ther- 
modynamic consistency. This means that in the limit of 
infinitely sharp interfaces between the different phases, 
the phase field equation of motion reduces to the phe- 
nomenological sharp interface equations that are adopted 
in a variety of analytical and numerical studies of mov- 
ing boundary problems. Demonstrating thermodynamic 
consistency for a class of systems places the usage of a 
phase field for this class upon a solid theoretical foun- 
dation [3|. Once thermodynamic consistency has been 
proven, a phase field model can be formulated for a prob- 



lem, for which sharp interface models are not yet avail- 
able. Bi and Sekerka [19| showed thermodynamic con- 
sistency for CH and AC systems. Simultaneously, they 
derived the equation of motion for such systems from an 
entropy functional. This is an alternative to applying 
Eqs.[3]or[3]to the free energy functionals[T]or[21 and this 
alternative has the advantage of being more generally 
applicable. 

So far, our discussion involved only systems with local 
equilibrium, where the movement of interfaces is 'slow', 
and classic irreversible thermodynamics applies. The 
movement of interfaces can be much faster in processes as 
e.g. rapid solidification 20|,|2l|, where the system can be 
far away from local thermodynamic equilibrium. Galenko 
and Jou [22| have shown how extended irreversible ther- 
modynamics can be used to develop a thermodynami- 
cally consistent phase field model that can deal with such 
fast dynamics using the entropy functional formulation 
(MCH and MAC models in Tab.]!]). Fast dynamics could 
be included by using an extended space of independent 
variables formed by the union of the classical 'slow' phase 
field variables, whose behavior is governed by conserva- 
tion laws, or an evolution of the phase field, and the 
additional space of 'fast' non-conserved thermodynamic 
variables, which consists of nonconserved variables with 
a relatively high rate of decay, for example, the fiux of a 
conserved phase field variable, or the time derivative of 
a non-conserved phase field variable [22l |. 

For systems with periodic stable states, Ste- 
fanovic et al. [18| were recently able to formulate an 
equation of motion that does not only describe time evo- 
lution at the 'slow' diffusive time scale - which is what 
the PFC equation of motion does - but also time evo- 
lution at the much faster time scale of elastic interac- 
tions. They called this the Modified Phase Field Crystal 
(MPFC) model, and demonstrated it provides physically 
meaningful results for nanocrystalling solids with high 
strain rates. However, its equation of motion was pos- 
tulated, and not obtained formally from a free energy or 
entropy functional. 

We will present now a single derivation that leads to 
the equations of motion for CH, MCH, PFC and MPFC 
models (sec Tab. |T| as special cases. Simultaneously, 
thermodynamic consistency will be proven for all four 
models. It is inspired both by the entropy functional 
formulation of Bi and Sekerka [19[ and the extended set 
of slow and fast thermodynamical variables by Galenko 
and Jou |22| . One distinction from previous unifica- 
tion attempts is that our approach does not require ad 
hoc assumptions of including the squared gradient and 
the squared Laplacian terms in the free energy func- 
tional pj, 12, [iSl, [iJI or in the entropy functional ISJ, \2^ . 

The derivation applies the extended space of indepen- 
dent variables E. We restrict our analysis to a conserved 
phase field variable tp and use the following extended 



space of independent variables E: 

E = {C} U {F} = {i'} U {J^,Vi^,V2^}, 



(5) 



Where {C} is the classical 'slow' phase field variables, 
and {F} is the extended space of 'fast' non-conserved 
thermodynamic variables. The flux j.^ was introduced in 
Ref. 122, and the two new extended independent variables 
Wijj and V^V in Eq. ^ are, to our knowledge, intro- 
duced here for the first time. With these two extended 
independent variables, it becomes a natural result to in- 
clude the squared gradient term [Vi/'P and the squared 
Laplacian term (V^'0)^ in the entropy functional, and 
obtain all four equations in a unified model. 

For the local nonequilibrium system described above, 
we postulate the existence of a local generalized entropy 
density s whose set of variables is the extended space 
E of Eq. ^. The generalized Gibbs equation for s is 
described by 

Siljjj^, VV, VV) = SeW + S„e(j*V: VV, VV), (6) 

where Sg is a local equilibrium contribution defined on 
the set {C} of the classic 'slow' variables, and s„e is a 
nonequilibrium part of the generalized entropy defined 
on the space {F} of the independent 'fast' variables, and 
can be written in the form 
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ai, a2 and a^ are system-specific parameters, and could 
be either positive or negative. The first term —ai\j,p\'^ /2 
is from Ref. 122. It has the physical meaning of a dif- 
fusion constant, and thus can be determined quantita- 
tively. The second term — a2|V'0p/2 and third term 
— q;3(V^7/')^/2 are obtained by symmetry arguments: 
since we treat all the additional extended variables j^, 
Vi/j, V^V on equal footing, there is no reason not to 
write all the terms in a similar way. The physical inter- 
pretation of the parameters a2 and a^ was discussed in 
Ref. llJ, and we will come back to this point later in the 
paper. 

For the system described by the extended set E of vari- 
ables, Eq. (O, the local balance law for the phase field 
variable ip is given by 



dtp f 

and the evolution of entropy density is defined by 



ds ^ ^ 



■ <^s, 



(8) 



(9) 



where as is the non-negative local entropy production. 
The entropy functional takes the following form: 



According to Eq. ([9]), the change of the total entropy S 
in time t is described by 



dS_ _ ldS_ 
m ^ \~dt 
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where 



dS 
dt 
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is the external exchange of entropy due to entropy flux 
is and 
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is the internal production of entropy due to dissipation 
within the system. In Eqs. (fT2|) and (fT3|) . A is the outer 
surface of sub-volume V and n is the normal vector to 
the surface. On the other hand. 
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and 



VV' • VV' = V ■ {tpVtP + V^tpj^) - v(vV) • U- 
Moreover 

v^v^v^v- = V • [v'^tpytp - tpv{y^tp) ~ vVJv] 
+v(vV)-j^. 

According to Eqs. ([TTl), dH]), and dni), we obtain 

ds ^ 
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+ Q3[V2^V^-V'V(VV)-VVJ0] (14) 
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Because as is non-negative, we get 
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Substituting j^ in Eq. ([8|) wc obtain 




aiM^ + ^ = V ■ MV ("aaVV - a^V^ - 


-!)<"' 


Recalling that 




ds 1 df 
dip Tdip' 


(18) 



to find models for a broader range of physical systems, 
for example, the yet-to-be-developed Modified Swift Ho- 
henberg (MSH) model in Tab. U as the analysis for a 
non-conserved phase field variable can be done in a sim- 
ilar way, with flux being replaced by time derivative in 
the extended space Eq. Q. 



and assuming that the mobility M is a constant, we have 
a.M^ + ^ = -V^ [a^TVH^ - a.TVH^ + ^ j(19) 



Eq. ^91) is the central result of the paper. It is now 
straightforward to sec that Eq. P^ is identical to the 
Modified Phase Field Crystal (MPFC) equation [3. If 
ai ~ 0, this equation reduces to the Phase Field Crys- 
tal (PFC) equation. If a^ = 0, it reduces to the Mod- 
ified Cahn-HiUiard (MCH) equation ^. Moreover, if 
ai = as = 0, it reduces to the Cahn-HiUiard (CH) equa- 
tion. Since «!, a2 and a^ are associated with the ex- 
tended independent variables j^, Vip and V^Vj we see 
the importance of including those variables in Eq. ([5]) 
and how each variable plays a different role in obtaining 
the existing models. 

By comparing Eq. \W\ with Eq. 25 in Ref. [23, wc find 
ai = Tu/M, where tjj is the relaxation time for the so- 
lute, td can be considered as the time needed to smooth 
concentration inhomogeneities by diffusion. As to the 
other two parameters a2 and as, a comparison with Eq.[T] 
and Eq.[2]revcals that a2 (related to the interfacial width) 
is positive when a^ ~ 0. When as ^ 0, a2 has to be neg- 
ative and as has to be positive to form a periodic phase 

In summary, we have shown how the four models from 
Tab. HI that have a conserved phase field variable, can be 
derived in a unified way by introducing two new extended 
independent variables in Eq. ([5]). This puts the MPFC 
model - which was originally based on a postulated equa- 
tion of motion - on a more solid theoretical foundation. 
Our derivation implies a proof of thermodynamic consis- 
tency for all four models in a unified framework. Mod- 
els with and without fast dynamics were treated on the 
same footing. As an outlook to further work, we suggest 
exploiting further the concept of thermodynamic consis- 
tency in the PFC case: derive sharp interface equations 
that correspond to the amplitude equations derived for 
PFC models [23|, |2J] for physically relevant problems, 
where a microscopic sharp interface model does not yet 
exist. As sharp interface models offer different possibili- 
ties for analytical treatment, they can contribute further 
to a more comprehensive understanding of the underly- 
ing physical system. In addition, the methodology pre- 
sented in the paper could be potentially further applied 
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